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In this work the Schmidt number of the two-photon state generated by parametric-down conver- 
sion (PDC) is evaluated in the framework of a fully spatio-temporal model for PDC. A comparison 
with the results obtained in either purely spatial or purely temporal models shows that the degree of 
entanglement of the PDC state cannot be trivially reduced to the product of the Schmidt numbers 
obtained in models with lower dimensionality, unless the detected bandwidth is very narrow. This 
result is a consequence of the non-factorability of the state in the spatial and temporal degrees of 
freedoms of twin photons. In the limit of a broad pump beam, we provide a geometrical interpreta- 
tion of the Schmidt number, as the ratio between the volume of the phase matching region and of 
a correlation volume. 

Introduction 

The process of parametric down-conversion (PDC) occurring in a nonlinear crystal is a widely employed source 
of entangled photons, which are ubiquitous ingredients in modern quantum technologies. An appealing aspect of 
this source is the possibility of generating high- dimensional entanglement, both in the sense that entanglement is 
generated in various degrees of freedom of the photon pair (polarization, time-energy, position-momentum), and 
because spatial and temporal entanglement is realized in a high-dimensional Hilbert space, due to the naturally 
ultrabroad bandwidths of the process. High dimensional entanglement is attracting more and more attention, because 
of its potential to increase the capacity and the security of quantum communication channels, and the precision of 
quantum metrological techniques P]. 

In this context an obvious question concerns the effective dimensionality of the entanglement of the PDC state (or 
alternatively the number of entangled modes generated by the process) , usually quantified by the so-called Schmidt 
number [H |3] . Traditional approaches typically concentrate on a single degree of freedom at a time, depending on the 
application considered. For example, the dimensionality of the temporal entanglement has been evaluated in various 
configurations, including spontaneous PDC (4»i5j, quantum frequency combs generated by a syncronously pumped 
optical parametric oscillator [SI [7], waveguided PDC [5]. The degree of transverse spatial entanglement of PDC 
[21 1101 115j is of paramount importantance for assessing both the dimensionality of the orbital angular momentum 
entanglement (see e.g. jll - 16 ) and the resolution of quantum imaging techniques |17| . In these studies a net 
separation between the spatial and temporal degrees of feedom was often justified by the assumption of a narrow 
filtering in the neglected degree of freedom. 

. However, as for many nonlinear optical processes, PDC is ruled by phase matching, which imposes an angular 
dispersion relation linking the frequencies and the angles of emission of the generated photons in a non-factorable way. 
This implies a strong coupling between temporal and spatial degrees of freedom, which has been recently evidenciated 
by the so called X- entanglement [T51I22j . a feature shown for example by type I bipotons in conditions close to collinear 
phase matching, whose temporal delays and the transverse spatial displacements at the crystal exit face are linked by a 
proportionality relation (corresponding to an X-shape of the spatio-temporal correlation in any plane containing time 
and one transverse coordinate). The space-time coupling is often regarded as a negative feature, because it affects e.g. 
the purity of the purely spatial entanglement when temporal degrees of freedom are neglected [23| . However, it also 
represents a valuable resource for engineering the quantum state of biphotons, since the spatial degrees of freedom can 
be used to tailor the temporal entanglement [181 119) in order e.g. to realize an ultrabroadband temporally entangled 
state [21]. 

In this work we adopt a fully spatio-temporal model of PDC, in order to investigate the effect of the non-factorability 
of the state in space and time on the dimensionality of the biphoton entanglement. We shall concentrate on the 
evaluation of the spatio-temporal Schmidt number of PDC entanglement in the ultra-low gain regime of PDC ^ 



^ The Schmidt decomposition of the state in the full spatio temporal domain is very hard to achieve (even numerically), due to the 
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Based on a general formula for the Schmidt number, which involves integrals in 12 and 6 dimensions, we both obtain 
numerical evaluations and analytical results, the latter being valid when the profile of the pump driving the process 
is broad enough. In the same limit, we introduce a useful geometric interpretation of our results, which shows 
that the Schmidt number quantifying entanglement is basically the ratio between the volume of the region where 
phase matching efficiently occurs and a correlation volume, thus being proportional to the number of spatio-temporal 
correlated modes. 

The same methods of analysis are used to calculate the Schmidt number in a purely 2D spatial model and in a 
purely ID temporal model of PDC. An important result that will be demonstrated is that the Schmidt dimensionality 
of the 3D model of PDC cannot be trivially reduced to the product of the Schmidt numbers in the models of lower 
dimensionality. As a matter of fact, the Schmidt dimensionality of the full PDC state is by far larger than what 
would be expected from an approach that considers space and time as separable degrees of freedoms, showing that 
an enormous number of entangled modes is available for the down-converted light. This result is again a consequence 
of the intrinsic non-factorability of the state in its spatial and temporal degrees of freedom, and shows that using a 
full spatio-temporal model for describing PDC is essential in order to correctly quantify the degree of entanglement 
of the state. 

The results presented in this work extend and complement those presented in [25] . where a 2D spatio-temporal 
model (1 temporal dimension -1-1 spatial dimension) for PDC has been investigated. 



I. STARTING POINT 



In [IHlIinim] the spatio-temporal quantum properties of the PDC light were described by studying the evolution of 
the quantum field operators throught the nonlinear crystal, and deriving input-output relations linking the operators 
at the crystal output face with those at the entrance face. Here we consider the equivalent state-formalism, in which 
the state evolves from the input to the output face of the crystal. We focus on type I PDC, in the regime of ultra-low 
gain where the probability of generating a photon pair in each spatio-temporal mode is small (more precisely, the 
probability of generating more than one photon pair in each mode is negligible). The output biphoton quantum state 
can thus be written as a generic superposition of the vacuum state |0) and of a state with two photons generated in 
all possible spatio-temporal modes: 

I^PDC) = |0) + / dwi I dw2C{iJ3^,W2)AHw^)A^{w2m (1) 



where A is the quantum field operators for the down-converted field, and w indicates the full 3-D spatio-temporal 
Fourier coordinate with the short-hand notation 

w = (g, VL) , (2) 

where q is the transverse component of the photonic wave- vector q = (q^, qy), with respect to the mean propagation 
direction z of the pump field, and — uj — ujp/2 is the temporal frequency shift from the central frequency of the 
PDC emission. The coordinate in the direct transverse space-time domain will be denoted by 

i^{x,t), (3) 

where x = {x, y) is the spatial coordinate spanning the transverse plane at the crystal exit face, and t is time, with 
the convention 

w-^^q-^-nt. (4) 

C{wi,'W2) in Eq.Q is the probability amplitude of generating a photon pair in the spatio temporal modes Wi and 
W2 , and can be determined by exploiting the equivalence with the field formalism developed in [T51 [ini [H] . In these 
references the biphoton amplitude was calculated in terms of the field-field correlation at the crystal output face 

^{W1,W2) = {A{wi,QA{w2jc)) (5) 

- Ap{wi -t- W2)smc ^ ' e' 5 (6) 



(27r) 
where 



hyperbolic geometry of the phase matching relations. However, the Schmidt number can be computed even without performing the 
decomposition. 
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- A^w, Ic) is the output field operator, and the expectation in Eq.([5| is taken on the input vacuum state; 

-g is the dimensionless gain parameter, proportional to the pump peak amplitude, the crystal length and the nonliner 
susceptibility; 



M^)--^ I lA^MOe-'^"' ^ (7) 



-Ap is the Fourier transform of the pump beam profile at the crystal exit face: 

(27r)3/2 

where normalization is such that Ap{S, = 0) = 1; 

-A is the phase matching function, which accounts for the conservation of longitudinal momentum in the microscopic 
PDC process 

A(wi, W2) = [ksziwi) + ksz{w2) - kp:,{wi + W2)] Ic (8) 

ksz being the longitudinal component of the (ordinary) signal wave vector, kp^ the analogous quantity for the 
(extraordinary) pump, and Ic the crystal length. 

We remark that the right hand side of Eq. Q is the first order term in the parameter g of a perturbative expansion of 
the full solution of the propagation equation of field operators in the nonlinear crystal, so that expression Q is valid 
only in the very low gain regime g <^ 1. Similarly, as it is well-known, the right hand side of Eq.Q shows the zeroth 
and first order terms in g of a perturbative expansion of the full PDC state. 

By using the equivalence between the two formalisms, and by calculating the field correlation on the generic output 
state ([1]) we also obtain: 

ip{wi,w2) = {il-'PDc\A{wi,0)A{w2,0)\tppr,c) 

= C{wi,W2) + C{w2,Wi) = 2C{wi,W2) ■ (9) 

where we used the symmetry properties of the state. Thus the two-photon state has the well known form 

IV'PDC) = |0) + ^/ / dW2^P{wi,W2)A\wi)A^ {W2)\0) (10) 

with the biphoton amplitude ■0 being given by ([6]). 

Apart from the biphoton amplitude, the other quantity of interest is the coherence function of the signal field, 
which after long but simple calculations can be derived from ([lO| as: 



G{w,w') :- {^pj,c\A\w)A{w'Mpj,c) (H) 

= / dw2'll^*{w,W2)i'{w' ,W2) (12) 



From this equation, the total number of PDC photons is obtained as 

N ^ dwi(V'PDc|^^(wi)^(wi)|V'PDc) 



dwi J dw2\ipi'Wi,W2)\^ (13) 

In type I PDC a signal and an idler fields are in principle not distinguishable. A bipartition of the system can 
be introduced in various ways, for example by sorting photons depending on their propagation directions (positive 
or negative with respect to any transverse axis) [22 , which, however, is a good bipartion only in the case of a very 
broad pump waist, when twin photons are always created with symmetric tranverse wave-vectors ±q. In the ultra-low 
gain regime where at most a single photon pair at a time is detected, a very general bipartition is obtained by simply 
considering the two output modes of a symmetric beam-splitter 

Ai{w) = ■^[A{w)+ia„{w)] 

Mw) = -^[iA{w)+ay{w)] (14) 
v2 
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where ay is a vacuum field operator. By substituting in the state ( 10 ) the inverse of ( 14 ) one gets two terms that 



describe creation of two photons into each of the output modes 1 and 2 of the beam-splitter, and a term that creates 
one photon in mode 1 and one photon in mode 2. Following the literature treating the degree of entanglement in 
PDC [H [51 [HI [TO] , we shall consider rather than the full PDC, the state vector conditioned to the measurement of a 
photon pair (the vacuum and two-photon terms are dropped) . Assuming that detectors are placed at the two output 
modes, and coincidences are detected, the state conditioned to the appearence of a coincidence takes the form (a part 
from global normalization factors) 



dwi I dv32Hwi,w2)A\{vSi)Al{w2)\Q)i\0)2 



(15) 



The degree of entenglement of such a conditional state has been investigated in previous literature in the purely 
temporal [H [S] or purely spatial [IHl [E] domains. 



II. THE SCHMIDT NUMBER OF PDC ENTANGLEMENT: INTEGRAL FORMULA 



A good quantifier of the degree of entanglement for continuous variable pure states is the so-called Schmidt number, 
defined as the inverse of the purity of the state of each separate subsystem 



K = 



1 

TV{pf} 



(16) 



where pi is the reduced density matrix of the subsistem 1. In connection with the Schmidt decomposition of the PDC 
conditional state, the Schmidt number is recognized to give an estimate of the number of Schmidt modes participating 
in the entanglement, i.e. of the effective dimensionality of the entanglement [10'. 



We will derive an integral formula for the Schmidt number in the case of the conditional state ( 15 ), similar to what 
obtained in Refs.[S], [TO]- First of all the state (151 is not normalized, 



= / dwi / dw2\'4>{wi,W2)\^ — N 



From the system conditional density matrix 

the reduced density matrix of the subsystem 1 can be calculated (Appendix |A| as 

Pi = Tr2{p} 



(17) 



(18) 



(19) 



Notice that in the limit where the coherence function becomes a Dirac-delta - i.e in the limit of a monochromatic 
plane-wave pump- the reduced density matrix becomes a sum of projectors onto 1-photon states. 
Next we calculate the purity of such a reduced state: 



1 

7V2 



dwi I dw'^\G{wi,w'i)\'' 



An integral formula for the Schmidt number can be threfore written as: 



K. 



B 



where 



B = 



dwi / dw'^lGlwijw'i)]^ 



dwi J dw2 J dw'i J dw2[i^{wi,W2)i^{w'i,W2) 

ij*{wi,W2)ij*(w[,W2)] 



(20) 



(21) 



(22) 



(23) 



and N is given by Eq.(17). 
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III. THE NEARLY PLANE-WAVE PUMP APPROXIMATION 



In order to evaluate the Schmidt number of the two-photon state from formula (21) we face the problem of calcu- 
lating the 6-dimensional and 12-dimensional integrals involved in the calculation of N and respectively . These 
integrations can be numerically performed, but in the following we will provide more trasparent results based on the 
approximation of a broad enough pump profile. 

Let us come back to the expression (|6]) for the biphoton amplitude that we rewrite as 

'lp'{wi,W2) = Ap{wi + W2)V{wi,W2) (24) 

N • A(wi,zZ;2) j ^("i."2) 
V{wi,W2) = smc^^ ^e* ^ (25) 

where A is the phase matching function defined in Eq. and we got rid of the costant g/ (27r)'^/^ that factors out in 
the ratio K. = N"^ / B. We now introduce the pump spectral coordinates Wp = Wi + W2 '■= {qp, i^p)- Provided that Up 
is the transverse waist of the pump beam at the output crystal face, and Tp its duration, the pump Fourier tranform 
Ap dies out on the scales 6qp — 2/ap, 6flp — 2/Tp. This claim is exactly true for a Gaussian pump profile 

Ap{x,t) = e-"'/"'e-*'/"p' (26) 

M<1P^^P) = "^e'^>'/'e'^nri/^ (27) 

The function V is strongly peaked along the curve where phase matching takes place. As elaborated in detail in 
the Appendix B of |19| . for a broad enough pump the variation of this function with respect to the pump spectral 
coordinates can be neglected. In other words, V{wi, —wi + Wp) does not vary significantly with Wp on the scale over 
which the pump Fourier profile dies out: 

Apiwp)V{wi,-wi + Wp) w Ap{wp)V{wi,-wi) 

Ap{wp)Viwi) (28) 

We call this approximation nearly plane wave pump (NPWP) approximation. It is based on making a Taylor expansion 
of V in power series of the pump variables Wp, and on finding the conditions under which the first order terms of the 
expansion can be neglected with respect to the zeroth order term [TH]. These conditions can be summarized as: 

Tp > TGVM = \—-—\ (29) 
Vgs Vgp 

dk 

C^p > Iwalk-off = l^c^^l- (30) 

Here tqym is the maximum delay time between the signal and the pump wave in crossing the nonlinear crystal, due 
to the mismatch between the group velocities Vgs, Vgp of the ordinary signal and extraordinary pump. Iwaik-off is the 
maximum lateral walk-off between the two waves, associated to tilt of the Poynting vectors In the example of a 4mm 
BBO crystal, pumped at a wavelenght Ap = 527nm, we have tgvm ~ 500fs, l^aik-off ~ 250^m, so that the NPWPA 
is within the reach of practical experimental generation of PDC photon pairs. 



Appendix [B] the integral formula (21) takes the form 

2 



The use of this limit simplifies remarkably the expression (21 ) of the Schmidt number. As reported in detail in the 

le form 

IdQApiipY'^ 



K. = 



[Jdw\V{w)\^Y ^g^^ 



B JdQApiipW i^^y Idw\V{w)\^- 
The integrals involving the pump coordinates can now be readily performed. By assuming a Gaussian pump profile 



as in ( 26 ) , we easily obtain: 



_ Jdw\V{w)\^ Jdw sinc^^ 
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As we shall see in the following, under rather general conditions ^ the first term at r.h.s. of Eq.(32) is a purely 
numerical factor, namely 

Jdw\V{w)\' ^ 3 

Jdw\V{w)\'' ^ 2' ^ ' 



Thanks to this circumstance, the result of Eq.(|32|) has a transparent geometrical interpretation. The term 
f dw sinc'^ A(^M)) interpreted as the volume of the region in the {q,fl) 3D space where phase matching oc- 

curs, since the sinc^ function has a sharp maximum where A((f, O) = (see Figjl]). This corresponds to the portion 
of the (q, fl) space where the probability of photon-pair production is not negligible. 

On the other side the quantity at denominator is (a part from numerical factors) the spectral volume of the pump 
SqpSilp = This quantity defines the correlation volume, i.e. the size the of the correlated modes, because 



the expression (24) tells us that in the NPWPA, the width of the biphoton correlation as a function of wi + 1^2 is de- 
termined by the pump Fourier profile. Thus Jq^^Op represents the uncertainty in the determination of the transverse 
wave vector 5*2 and frequency ^2 of a photon, once the transverse wave vector qi and frequency Vti of its twin have 
been determined. Therefore K. is on the order of the ratio: 

volume of the phase matching region 

IC oc , — -. , (34) 

correlation volume 



and can be interpreted as the number of correlated modes partecipationg to the state. The formula (34) gives us a 



simple geometric interpretation of the Schmidt number, which will be very useful in understanding some results. 




FIG. 1: Geometrical interpretation of the Schmidt number. The figure plots a) the phase matching region in the {qx, f2) plane, 
here defined by |y((f, > 0.1 (the full 3D volume has a biconical shape), and b) the correlation volumes, which in the 3D 
picture would be Gaussian bullets of size determined by the spectral estension 5q^5Q,p of the pump. Case of a 4 mm type I 
BBO crystal, pumped at Ao = 527nm for coUinear phase matching. 



^ this conditions amount to requiring that phase matching occurs within the spectral region considered, a conterexample being provided 
in SectiorfVll 
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IV. SCHMIDT NUMBER OF 3D X-ENTANGLED BIPHOTONS 



We can proceed further, and find an analytical approximation of the Schmidt number result of Eq.( |32[ ) in the 
NPWPA. To this end, we need to calculate integrals over the phase matching curves of the form J dw\V {'w)\'^ , 
J dw\V{w)\'^. Our calculations are based on the use of two further approximations (in additions to the NPWPA): 
i) The first approximation consists in a quadratic expansion of the phase mismatch function with respect to q and fJ, 
equivalent to adopting the paraxial and quadratic dispersion approximations: 



A{q,n) ^ Aolc - 



fi2 
S2q 



(35) 



where Aq — 2ks — kp is the coUinear phase mismatch at degeneracy, and 



no 




kglc 



(36) 
(37) 



with kg = fcs(O), fc" = cPks/(m?\o. This expansion is strictly valid only for small ft (close to degeneracy) and small q. 
For the remaining of this section we assume conditions of coUinear phase matching, AqZc ~ 0, where the phase matching 
curve in the plane {q,fl) has the characteristic hyperbolic geometry shown in Figjl] The parameters qo, associated 
to spatial diffraction and Oq, associated to group velocity dispersion (GVD), define the typical variation scales of 
|y(w)P along q and fi, respectively. In the example of the 4 mm BBO crystal, their values are go ~ 5 x 
ilo w 0, 76 X 10"Hz. 

ii) The second approximation consists in substituting the sinc^ ^ A(t») ^ with a box function, with the same value of 
the indefinite integral 



smc 



e(- 

elsewhere 



(38) 



which satisfies J sinc^{x)dx = f Xa{x)dx = n. Here the parameter a can be used in principle as a fitting parameter. 
Approximation ( |38[ ) seems very rough, but it turned out surprisingly accurate: the rationale behind this result is that 
the sinc^ has a sharp peak where A{w) = 0, and in order to evaluate its integral in the 3D space it is more important 
to take into account the geometrical shape of the curve where its maximum lies rather than the detailed shape of the 
peak. Figure [2] compares the box function approximation to the true phase matching function in the example of the 
BBO crystal. Here, substantial deviations appear at large values of and q because of the failure of the quadratic 
approximation for phase matching. 

Notice that if we also assume that |y(?i;)|^ can be approximated by the box function XqI^^)? with the request that, 
this time. 



J dxXa{x) = — = J dx sinc^(a;) 



we find the correct value of a 



Itt. This also shows that, within these approximations, the ratio 



/ dw\V{w)\^ _ / dwxa{w) a 3 
Jdw\V{w)\^ Jdwxliw) ~ ^ ^ 2 



(39) 



(40) 



The box function approximation allow us to evaluate easily the integrals over the phase matching curves inside Eq. 
( 32 1 . In this evaluation, we assume that our model describes a measurement performed over a large but limited 

). For simplicity, here we do not pose limits to the spatial bandwidth (which 
After some calculations we obtain the two different results depending 



spectral bandwidth fl e (—flrnax,^ 
will be instead done in the following section) 
on the detected bandwidth Qm„T = ^n°^ : 



• small bandwidth result {flmax < y/ct) 



2n 2 



(41) 
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box ftiactiou 
a ~ 1.5tz 




1 2° 2 . 




FIG. 2: a) Box function approximation (381 of tiie sinc^ function, b) Comparison between the true phase matching function 
r2)P and its box function approximation (boundaries of the box function are shown by dashed lines). Cohinear phase 
matching case (AqIc = 0, 6p = 22.934°), Ic = 4mm. 



• large bandwidth result {ftmax > \/«) 



IC = 



p 'p 



^7} 



la 



(42) 



First of all we observe that the condition on the bandwidth can be roughly reformulated as ^max being smaller or 
bigger than the characteristic GVD bandwidht f^o (since a is on the order of unity). The small bandwidth case 
corresponds to the situation where the portion of PDC emission intercepted by the measurement lies within the 
central region of the phase matching curve (see Fig[T]), where the phase matching has no hyperbolic structure. Eqs. 
(41 1, (42) tell us that in both cases the Schmidt number is proportional to the number of modes contained in a unit 



volume of the phase matching region: 



K, oc 



(43) 



However, as the detected bandwidth increases beyond the GVD bandwith Hq, the hyperbolic geometry of phase 
matching enters to play, and the Schmidt number result of Eq. ( 42 1 shows a linear increase with the bandwidth. 



We remind that the analytical expressions ( 41 1 and ( 42 ) estimate the Schmidt number within the NPWPA and 
the quadratic approximation, expressed by the condition ( 28 1 and ( 35 ) , respectively. In order to verify its validity. 



and at the same time provide a numerical estimation of JC in regions of the parameter space where the NPWPA does 
not hold, we performed a numerical evaluation of the general expression of JC given by Eq.(21),(13) and ([23|). As it 



involves a 6-dimensional integral for the evaluation of N and a 12-dimensional integral for the evaluation of B, the 
use of a Montecarlo integration is mandatory. We used the well known method of importance sampling |26j with the 
aim of improving the efficiency of the Montecarlo algorithm by increasing the density of the sampled points where the 
functions under the integrals are larger. A natural choice has been to sample some of the Fourier variables (namely 
the "pump" variables in Eq.(B5)) according to Gaussian distributions coincident with the Gaussian pump spectral 



amplitude ( 27 ) , which in the NPWPA represents the narrowest factor of the the biphoton amplitude in Eq. ( 24 1 . The 
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inplemented algorithm, very efficient for narrow specral pump profiles, allows the evaluation of B and N even in the 
region where the the NPWP approximation fails. No other substantial approximations are introduced, as the phase 
matching function is here evaluated by means of the empirical Sellmeier formulas [27) . 
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FIG. 3: Schmidt number results. Comparison between the analytic formula (41 1, (42 1 (solid red line) and Montecarlo simula- 
tions, without any approximation (square) and with quadratic approximation (triangles), a) Pump parameters are within the 
NPWPA b) Focused pump, beyond the NPWPA. CoUinear phase matching (Aq/c = 0, Op = 22.934°), Ic = 4mm. 



Fig. p^compares the analytic result (41 ), (42 ) with the Montecarlo numerics, performed both without approximation 
(squares) and with the quadratic approximation for phase matching (triangles). Part a) of the Figure is plotted for 
parameters of the pump within the NPWP approximation (althought very reasonable for an experimental realization) 
and shows an excellent agreement between the analytical curve and the numerics, in the range of validity ogf the 
quadratic approximation. Indeed, the analytic result follows very well a Montecarlo simulation performed with the 
quadratic approcximation, showing that the box function approximation captures the basic geometrical properties of 
the phase matching function. In plot b) the pump beam is more focused and as expected the analytic result deviates 
from numerics because of the failure of the NPWP approximation. 

Montecarlo calculations permit to obtain results also in the region of parameters beyond the NPWP approximation. 
Figure Rl plots the Schmidt number as a function of the parameter /3 — SqpSD,p/qQQo. The NPWPA approximation is 
valid omy for /3 <C 1, i.e. when the widths 6flp, Sqp of the pump Fourier profile are much smaller than the characteristic 
scales of variation of phase matching rig, 9o, respectively The Montecarlo result shows a decrease of the Schmidt 
number as 1//3 for /3 <C 1, as predicted by the analytic result (42) in the NPWPA (red solid line in the Figure). 
However, after reaching a minimum the Schmidt number increases again almost linearly with /3. This behaviour is 
very similar to the one predicted in a purely spatial model of PDC in Ref.|^ and can be understood as follows: for 



^ Actually the limits of validity of the NPWPA expressed by \29\ and ( |30[ l are typically much more restrictive than /9 <C 1- For 
example for the BBO crystal here considered Sqp < qo implies roughly ap > 2l^iff = 2/qo ft: 40/im, while <5f2p < SIq implies 
Tp > 2tqvd = ~ 26fs. 
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FIG. 4: Schmidt number as a function of /3 = SgpSfip/q^Qo- The blue squares plot the result of a Montecarlo calculation, 
without any approximation, and shows that IC after reaching a minimum increases again almost linearly woth f5. The red solid 
line is the analytic result, decreasing as 1//3 (see Eq.(42|, valid only within the NPWPA (small /3). CoUinear phase matching 
(Ao^c = 0, 6lp = 22.934°), Ic = 4mm. 



a broad pump, when the NPWPA is valid, the width of the correlation is determined by the pump Fourier profile, 



and the number of spatio temporal modes can be estimated as in formula (34) as being proportional to the volume 
of phase matching divided by the pump spectral volume K, oc q^nQ/SgpSflp = For a very focused pump, instead, 
the phase matching function in the q direction has a smaller scale of variation than the pump Fourier spatial profile, 
so that the width of correlation is rather determined by the characteristic width go of phase matching, and we can 
suppose that the number of modes is now /C cx Sq^/q^ (x /3. 



FACTORABILITY OF THE SCHMIDT NUMBER IN ITS TEMPORAL AND SPATIAL 

COMPONENTS 



In the literature the Schmidt dimensionality of twin photons is often calculated within models of PDC restricted to 
the spatial or the temporal domain (see (H El O [10] ■ The non-factorable character of the spatio-temporal correlation 
demonstrated in [18l [19] suggests us that the full 3D spatio-temporal Schmidt number is not trivially given by the 
product of the spatial 2D and the temporal ID Scmhidt numbers. In this section we would like to understand this 
point. 

To this end we consider models for PDC in lower dimensionalities, and follow the same procedure outlined in 
the previous sections to calculate the Schmidt number. The purely spatial 2D model is obtained by neglecting the 
temporal coordinate and setting f2 = 0. Similarly the purely temporal ID model neglects the spatial coordinates and 
sets q = Q. The starting point of the analysis is in both cases the general integral formula for the Schmidt number 



(21 1 where we have now to interpret the Fourier coordinates as: 



n en in id 

w = < q en"^ in 2D 
q,n en^ in 3D 

Simlarly, in the expression involving the coordinates in the direct space: 

t en in ID 
^ = { X en"^ in 2D 
x,t en^ in 3D 



(44) 



(45) 



For example, by introducing the NPWP approximation in the various models ( clearly NWPA in the spatial model 
means that the pump has a broad waist, while the temporal model assumes a long enough pulse duration), we obtain 
the NPWPA expression for the Schmidt number in an arbitrary D-dimensional model: 



/C = 



[J dw\V{iv)\^]^ 
i2n)D Jdw\V{w)\^ 



(46) 
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A. Spatial Schmidt number IC2D 

By performing calculations similar to those reported for the 3D model, we derive an expression for the Schmidt 
number in the purely spatial case, valid within the NPWPA and the quadratic approximation for phase matching. 
The latter one corresponds to approximating the phase matching function as 

„2 



^2D{q) = ks{q) + ks{-q) - kp 



In the 2D case the result depends on the spatial bandwidht q^^^ = Qmax/Qa intercepted by the measurement: 



(47) 

(48) 
(49) 

where we remind that a ss I.Stt. In Figj5]this curve is compared with an exact Montecarlo calculation performed 



1^2D = 



< 



> 




-■- K 2D Montecarlo 

K aiial}lic 

a= 1 . 5 71 a =600 iim 
p 



6 8 10 12 14 



FIG. 5: 2D spatial Schmidt number IC2D as a function of the collected spatial bandwidth qmax- The solid red line is the analytic 
result (48 1, (491, the squares plot the Monte Carlo numeric result. The waist of the pump beam is Op — 600/im. 



in the 2D model. Beyond noticing that the two results agree qualitatively, we remark that, differently from the 3D 
case, the 2D Schmidt number saturates to the maximum value JC2Dmax = I^^'b^o ~ ^''^lol^^^y This behaviour can 



be explained with the help of the geometrical interpretation (34), valid in the NPWPA, which evaluates the Schmidt 



number as the ratio between the volume of the phase matching region and the correlation volume. In the 2D case 
phase matching is described by Eq. (47), so that in the {qx,qy) plane phase matching occurs within a circle of area 
~ TT^. For increasing qmax, the PDC photons are detected in increasing large circular regions, so that the Schmidt 
number increases quadratically with qmax until the border of the phase matching region qmax = qoV^ is reached. 



B. Temporal Schmidt number /Ciu 

We now consider the purely temporal model of PDC, by setting q = 0. We perform again analytic calculation in 
the NPWPA and the quadratic approximation for phase matching, that in the ID case reads 

Aiom ^ ks{n) + ks{-n) - kp ^ (50) 

"0 

The analytic expression for JCm in these limits, obtained by using the box function approximation, depends on the 
collected temporal bandwidth ilmax = ^max/^a'- 

IClD = \^Tp^O^^^ ^max < (51) 
K-iD = ^Tpilo i^max > \foL (52) 

Figure |6] plots this analytical result toghether with the exact Montecarlo ID calculation. Also in this case, similarly 
to the 2D case, the Schmidt number saturates to the maximum value ICiumax ~ TpQQ, because phase matching occurs 
only inside an interval of size « Jloi so that by increasing ^max beyond the critical value -^/a^o j '^6 begin to consider 
regions where there is no phase matching, which do not contribute to the integral. 
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FIG. 6: ID temporal Schmidt number K\d as a function of the maximum temporal frequency ^max- The soUd line shows 
the analytic result (within the NPWPA and quadratic approximation), the squares plot the numeric exact result. Pump time 
duration Tp = Ips. 
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FIG. 7: Gomparison between the Schmidt number K, calculated in the full spatio-temporal model and the product of the 
Schmidt numbers YZ,\d x K,2d, calculated in the purely ID temporal and 2D spatial models. The abscissa is the collected 
temporal bandwidth, set equal to the spatial bandwidth in normalized units. Lines plot analytic results, symbols provide the 
Montecarlo results. 



C. Comparison 

We are now in conditions of comparing the results obtained in the models of various dimensionalities. To this end. 



we have to reformulate slightly the 3D problem. In Sectior IV we calculated /C as a function of the collected temporal 



bandwidth by assuming that no selection on the spatial bandwith was performed, i.e. q^nax = This is a possible 
correct choice to present results, but in order to have a clean comparison with the 2D and ID models, we need also to 
limit the detected spatial bandwidth. The simplest possibility is to set q^nax — ^max- This choice is justified by the 



fact that in the quadratic approximation, phase matching is realized along the lines — = zkS- so that by increasing 



imax 

go f^o 

simultaneously the spatial and temporal bandwidth 2^^^ = ^^ff^ one follows the phase matching curve. 
With this in mind, analytical calculations in the PWPA, quadratic approximation for phase matching and box function 
approximation can be performed. The result for the 3D Schmidt number is plotted in FigjT] (dashed line) toghether 
with the Montecarlo exact result (triangles) . In the same figure we compare these 3D results with the product of the 
Schmidt numbers obtained in the models with lower dimensionalities i.e. /C2rf x K^iD- From these plots it clearly 
emerges that the factorizability holds only when the detected bandwidth is small, i.e. when both qmax and ^Imax 
lie within the phase matching bandwiths qo, flo, respectively. However, as the detected bandwidth gets larger, the 
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FIG. 8; Colormap plot of \V{q,^)\^ in the non-collinear phase matching case, Aq/c > 0,(Aoic = 23.38). The dashed lines are 
the boundaries of the box function, where A(g*, fi) = ±q 



result in the fully 3-D spatiotemporal model grows linearly with the detected bandwith, and diverges clearly from the 
product K^d X which saturats to a fixed value oc q^^Q. This result can be easily understood with the help of the 
geometrical interpretation of the Schmidt number as the number of entangled modes contained in the phase matching 
region: close to the degeneracy and to the collinear emission, the phase matching region can be seen as a spherical 
region, which obviously factorizes in the spatial and temporal components, so that the number of spatio-temporal 
modes is trivially the product of the numbers of spatial and temporal modes times. Conversely, if the collected 
bandwith is large enough, than the biconical, non-factorizable geometry of phase matching comes into play, so that a 
full 3D model has to be used to correctly compute the number of spatio-temporal modes. 



VI. ENTANGLEMENT WITHOUT PHASE MATCHING 



Up to now we considered the case of collinear phase matching (Aq/c = 0), characterized by the fact that exact 
phase matching i^(w) = can be realized for any value of the transverse wave-vector of the photon pair, such that 
q/qQ — ±$7/^0. However, if the crystal is tuned away from the collinear conditions (Aq/c 0), there exist regions of 
the (g, 17) space where phase matching does not occur at all. In these regions, the probability of emission of photon 
pairs is low, although not zero. 

In particular, we focus on the case of non-collinear phase matching (Aq/c > 0). Figjs] shows the typical behaviour of 
the phase matching function sinc^(A(iZi)/2) . From the quadratic expansion of the phase mismatch 



A(g,f7)« Ao^c-^ + S' (53) 



we immediately realize that phase matching does not occur for modes close to collinear emission, namely having 

We have calculated the 3D spatio-temporal Schmidt number in these conditions: Fig. |9]plots the result as a function 
of the collected spatial bandwidth qmax- We find that for small bandwidths, such that the collected photons are not 
phase matched, the degree of entanglement is very high , actually larger than in the region where phase matching is 
realized. This result is apparently paradoxical, because in the absence of phase matching the probability of emission 
of a photon pair is very low and one could infer that the state should be very close to the separable vacuum state. For 
comparison Fig |10| displays the corresponding mean number of photons, which as expected, is indeed very low where 
there is no phase matching. 

However, in order to understand the results of Fig. [9] we have to remind that we are studying the degree of 



entanglement of the state (15), conditioned to the detection of a photon pair. This means that our calculation of the 
Schmidt number has lost track of the presence of a large vacuum contribution to the original PDC state, and the 
result has to be interpreted as photon pairs are emitted with very low probabilty, however, when a pair is detected it 
is highly entangled. 

The point to understand is therefore why the non-phase matched photon pairs are more entangled than the phase 
matched ones. The Schmidt number K. = ^/J^j-^j^ provides an estimate of the number of significant eigenvalues of 
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FIG. 9: Schmidt number K, in tlie non-collinear phase matching case , as a function of the collected spatial bandwith q, 
normalized to the diffraction bandwith go- Parameters are ap = 600pm, Tp = Ips 
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FIG. 10: Number of photons A'' in the non-collinear case, as a function of q^nax- Parameters are o-=600/^m, Tp = 1 ps, g = 0.001. 



the Schmidt decomposition, i.e. the number of entangled eigenmodes that participate to the modal decomposition. 
In the region where no phase-matching occurs (for q^^ax < -s/AqZc)) ^ can be in practice very large, since there 
is no mechanism for modal selection and all the modes in the collected band'width participate equally to the PDC 
process even though with a very low occupation probability. By contrast, when increasing q^ax towards -v/Aq/c and 
entering the phase-matching region, a strongly reduced number of phase-matched spatio-temporal modes contribute 
to the Schmidt decomposition, the few which are close to satisfy the phase-matching condition q = \/AoI^, 57 = 0. 
The number of significant eigenmodes is therefore reduced because phase matching operates a selection of the spatio- 
temporal modes that efficiently participate to the entanglement of the state. In other terms, the phase-matching 
mechanism favors only a small number of modes, reducing thereby drastically the dimensionality of the entangled 
state. By increasing q^ax above -v/Ao'c the Schmidt number K, starts again to increase, according to the geometrical 
interpretation of the Schmidt number as being proportional again to the phase matching volume 



VII. CONCLUSIONS 



In this work we have calculated the Schmidt dimensionality of the two-photon state generated by PDC in the 
ultra-low gain regime. Results have been produced with different degrees of approximation: Montecarlo calculations 
allow us to obtain the Schmidt number without relevant approximations, while when the pump beam is broad enough 
(NPWPA) we have demonstrated a trasparent geometrical interpretation of the Schmidt number, which can be seen 
as the number of entangled modes contained in the region where phase matching occurs. 

The same calculations have been performed in models restricted to the purely spatial or temporal degrees of freedom 
of biphotons. A remarkable result demonstrated is that, when the collected spatio-temporal bandwidth is large enough. 
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the Schmidt dimensionahty of the full spatio-temporal state cannot be trivially reduced to the product of the Schmidt 
numbers characterizing the entanglement in lower dimensions. Therefore, obtaining the Schmidt number in the full 
3D model is not a mere exercise of calculus: in order to correctly characterize the entanglement of twin photons it is 
necessary to consider simoultaneously their spatial and temporal degrees of freedom. This result is a clear consequence 
of the non-factorability of the state in apace and time, and mirrors the findings described in Refs. [TH], [H] where the 
spatio-temporal correlation of the biphoton state was shown to have a non factorable X-shaped geometry. 

The non-factorability has been demonstrated in this work only in the NPWPA. Question is still open whether in 
the opposite limit, that is, for a very focused and short pump pulse, the Schmidt number keeps the non-factorable 
character, and it is obviously linked to the more general question whether the state appears to be factorable in space 
and time in these conditions. We however remark that reaching this limit is in practice very demanding, because 
it requires that the pump Fourier profile at the exit face of the crystal is much broader than the width of the 
phase matching function. For a few millimeter crystal, this implies a pump waist on the order of tens of microns, 
and a pulse duration as short as few tens of femtosecond. While the first condition could be in principle achieved by 
strongly focusing at the end face of the crystal, the second is more demanding because of dispersion occuring inside 
the nonlinear material. 

An intriguing finding is that in the absence of phase matching, where the probability of emission of photon pair is 
very low, the Schmidt dimensionality of the state is huge, actually larger than in the regions where phase matching 
occurs. This counter-intuitive finding has been explained throught the mode selection mechanism performed by phase 
matching, which reduces the available number of spatio temporal modes. However, in evaluating this result and its 
usefulness for applications, one has to remember that the Schmidt number here analysed does not refer to the full 
PDC state, but to the state conditioned to the detection of a photon pair. 

Acknowledgments 

We acknowledge support from the Fet Open programme of the EC within the grant 221906 HIDEAS High Dimen- 
sional Entangled Systems 

Appendix A: Partial density matrix 



In order to calculate the partial trace of the density matrix (18) it is enough to consider the vacuum state [0)2 of 
the "idler" subsystem 2, plus the continuous set of 1 photon states 

{4{w2)\0),} (Al) 



It is convenient to write the density matrix ( 18 1 as 



p = |0) (0| M (A2) 
where M is the operator that annihilates a photon pair in any spatiotemporal mode (weighted by 4') 

A/ = y dwi J dW2'ip* {wi,W2)Ai{wi)A2{w2) (A3) 

We then have 

Pi = Tr2{p} 

= ^ { 2(01 10)2 |0)i i(0| 2(01 M 10)2 + J dw3 2(01 A2{w3)M^ 10)2 |0)i i(0| 2(0| MAUw^) \0)^'^ (A4) 
= — |y dw^ J dwi J dw2 J dw[ J dw2'ilj{wi,W2)i'* {w[,W2) 

X 2(01 A2(«J3)4(^^2) 10)2 4K)|0)i i{0\ Ai{w[) 2{0\A2{w2)AI{w3)\0) 2^ (A5) 

= dwi J dw2 J dw[i>{wi,w2)riw[,W2)A\{wi)\0)^ i{0\ Ai{w[)^ (A6) 

= ^{[dwi f dw[G{w[,Wi)A\{w,) |0), i(0| Ai(^)| (A7) 
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where in passing from ( A5 1 to ( A6 ) we used the relation 

dfi\A2{w2)A\{w'.^\Q).,=5{w2^lS'2) 

which conies directly from the commutation rules of bosonic operators. 



(A8) 



Appendix B: Derivation of formula (31) 



We wish here to simplify the general formula (21) by exploiting the NPWPA introduced in Section III (see Eq.(28)) 
We rewrite formula (21 1 as 



13' 



(Bl) 



with 



B' = j dwi J dw2 J dw'i J dw2Ap{wi + w2)Ap{wi + 'W2)Ap {W1+W2) 

Ap {w'i+W2)Ap {Wi + W2) ■ V{wi,W2)V{w[,W2)V*{w[,W^^ 
N' = I dWldW2\ApiWl+W2)\^\ViWl,W2)\'^■ 



m 

(B3) 



where we inserted the explicit expression (24) for the biphoton amplitude (without a costant coefficient). Here the 
function V depends on phase matching and is given by Eq.(25), while Ap is the Fourier profile of the pump. We start 



by simplifying the integral B' in Eq. ( B2 ) , by introducing the change of variables 

1 + W2,w'p w[ + 



Wi — Wi, Wp 



{Wi,w[w2,w'2) (wi,S 

With this change, B' becomes 

B' ~ J dwpdWpdSdwiAp{wp)Ap{wp)Ap {wp~S)Ap (w'p + &)V{'Wi,~wi + Wp) 
V(wi — 5, —wi + S + Wp)V*{wi — S, —wi + Wp)V*(wi, —wi + S + Wp 



(B4) 



(B5) 



In this expression the variables Wp and w' die on the scale of the inverse of the pump waist/duration, i.e. {2/ap, 2/Tp). 
Because of the presence of the terms Ap(w[ — Wi+ Wp) and .4* (wi — w[ + Wp), also the variable 6 — Wi — w[ dies out 
on the same scale. We can then make use of the NPWP approximation, which amounts to substituting 



V{wi,—wi+Wp) V{wi — S,—wi + S + Wp) 
sa V{wi, —wi) := V{wi) 



(B6) 



where we took into account that all the pump variables Wp, S, S + w'p dies out on the fast scale of the inverse of the 
pump waist/duration, over which the function V remains constant. Similarly 



V*{wi — + Wp) w V* {wi, -wi + 6 + Wp) 

« V*{wu-wi) -.^V^wi) 

Within the NPWP approximation we hence obtain a new expression for B' , which reads: 



(B7) 



B' ^ J dwpdWpd6Ap{wp)Ap{Wp)Ap {wp-S)Ap {Wp + S)x J dwi\V{wi)\'^. (B8) 
The integral over the pump variables can be further simplified by noting that it involves two convolutions: 

dwpAp{wp)A;{wp ± 6) = J dip\Ap{Q\\-^'i^--' = (27r)i [T{\Ap\^)] (±6), (B9) 

where the symbol J-{f) denotes the Fourier transform of the function f . By performing also the integration over 6 
we obtain 



d6[:F{\Ap\')] {5)[J^{\Ap\')] i-S) 



dip\Ap{ip)\\ 



(BIO) 
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~ 2 2 

where we used the Plancherel theorem / (Pw\f{w)\ — J <P^\f{£.)\ ■ This leads to: 



B' = (27r)3 



dw\Vpw{w)\'^ 



In order to complete the Schmidt number calculation we also have to evaluate N' , proportional to the mean 
number. With the usual change of variables {wi, W2) {wp = Wi + W2,wi), equation (|B3[) becomes: 



N' = / dwp \Ap{wp)\^ / dwi \V{wi,~wi+Wp)\' 



In the NPWP limit (see Eq. [28|) we get 

N' = J dwp\Ap{wp)\^ J dw\Viw)\^ , 
which, using the Plancherel theorem in the first integral, can be also written as 

N'= [ dQAp{Q\' I dw\V{w)\\ 



In NPWP limit the Schmidt number takes therefore the symplified form 

2 



K = 



~B 



I d^p\Ap{^p)\' 



[Jdw\Viw)\^^ 



Jd^ApiQl"^ (27Tr J dw\V{w)\^ 



(Bll) 
photon 

(B12) 
(B13) 
(B14) 

(B15) 



where the integrals now factorize into the pump and signal degrees of freedom. 
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